A quasilocal calculation of tidal heating 
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We present a method for computing the flux of energy through a closed surface containing a gravitating 
system. This method, which is based on the quasilocal formalism of Brown and York, is illustrated by two 
applications: a calculation of (i) the energy flux, via gravitational waves, through a surface near infinity and (ii) 
the tidal heating in the local asymptotic frame of a body interacting with an external tidal field. The second 
application represents the first use of the quasilocal formalism to study a non-stationary spacetime and shows 
how such methods can be used to study tidal effects in isolated gravitating systems. 



I. INTRODUCTION 

In many physical problems in gravitation, one is interested 
in the interaction of a nearly isolated gravitating system with 
an external universe. The interaction effects are computed in a 
"buffer zone" (see Sec. 20.6 of Ref. and Sec. IB of Ref. [|]) 
surrounding the gravitating system, in which the radius of cur- 
vature, scale of inhomogeneity, and rate of change of curva- 
ture are much smaller than the size of the body. The formalism 
of Thorne and Hartle [^] and Zhang has been used recently 
by Purdue and Favata to compute the gauge-invariant 
heating of a body interacting with an external tidal field. 

Until now, calculations of the sort described in these refer- 
ences have made use of pseudotensors to compute energy and 
momentum fluxes. However, quasilocal methods should be 
equally applicable in situations with a reasonably well defined 
buffer zone — in this case, the quasilocal surface can be con- 
veniently located in the buffer zone. While quasilocal meth- 
ods are not fundamentally different than pseudotensor meth- 
ods [^], an advantage of quasilocal method is that all quanti- 
ties (e.g., energy fluxes) can be computed in terms of real ten- 
sors on the quasilocal surface. Gauge ambiguities in the total 
amount of energy and energy flux such as those reported in 
Ref. and discussed in Ref. [Q] still exist for the quasilocal 
methods, but now the ambiguities can be understood in terms 
of distortions of the quasilocal surface and so their geometric 
origin is identified. 

In this paper we present a quasilocal formalism for com- 
puting the work done on a gravitating system by an ex- 
ternal universe. Our formalism is based on the quasilocal 
mass of Brown and York |Q| — the on-shell value of the grav- 
itational Hamiltonian — which coincides with the Amowitt- 
Deser-Misner energy at spatial infinity and the Trautman- 
Bondi-Sachs energy at null infinity ||7|||]. It is complementary 
to but independent of |^] which studied how motion of the ob- 
servers affects the Brown- York energy. We use our expression 
for energy flux to compute (i) the energy lost in gravitational 
radiation from a gravitational system and (ii) the heating of a 
body through interactions with an external tidal field. Problem 
(i) demonstrates that the formula for the work reproduces the 



known gravitational radiation flux formula when the quasilo- 
cal surface is located in the wave-zone. Problem (ii) repro- 
duces the calculation of Purdue using quasilocal methods 
and shows how these methods are applicable for problems in 
which the quasilocal surface is located in a buffer zone. 



II. QUASILOCAL ENERGY FLUX 

In this section, we derive an expression for the energy flux 
through a closed two-surface suiTounding a gravitating sys- 
tem. Our analysis closely follows Sec. V of Ref. which 
derives a conserved measure of mass for stationary systems. 
We relax the requirement that the quasilocal two-surface time 
evolution vector be a Killing vector of the spacetime and 
thereby obtain an expression for the rate of change in the mass 
of the system. 

Consider a gravitating system separated from the external 
universe by a (2 + l)-dimensional timelike boundary B. This 
boundary has an outward "radial" normal vector n", a metric 
7afc — dab — naTib induced by its embedding in the spacetime 
with metric gab, and an extrinsic curvature Qab = —^$nlab 
(with trace 8 = j^'^'Qab)- Let be the derivative opera- 
tor compatible with the metric •jab- Foliate the boundary B 
into closed two-surfaces ilt of constant time t; then the time 
evolution vector t°- on B satisfies f^Aat ~ 1 and can be de- 
composed into a lapse function N and a shift vector V"" on 
ilt via <° — Nu'^ + V"-, where is the timelike normal to 
fit embedded in B. The closed, spacelike, two-surface fit 
has an induced metric Oab = lab + UaUt and, viewed as a 
two-surface embedded in a three-dimensional spacelike hy- 
persurface S locally defined such that G TS, the extrinsic 
curvature of fit is kab = —^SnC^ab- A full discussion of the 
geometry of the boundaijy B and its foliation (including a dia- 
gram) may be found in [^. The notation there is substantially 
the same as here though w° is written as u". 

The Codazzi identity. 
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where = (67°'' — 8"^) / Svr, relates the extrinsic curvature 
of B to the spacetime Ricci curvature Rab- It then follows 
from the Einstein field equations that 



ab\ 
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We restrict our attention to a vacuum spacetime in which the 
stress-energy tensor Tab vanishes. Then, if <° is a Killing vec- 
tor field of the boundary metric ^ab, Eq- (H) is a conservation 
equation and the quantity 



(3) 



is a conserved measure of the total mass contained within 
the boundary It is the "non-orthogonal" Brown- York 
mass [^]^], up to a subtraction term that is required for it to be 
bounded for large surfaces in asymptotically flat spacetimes 
(see, e.g., Lau or Mann [|lo|]). 

When i° is not a Killing vector of the boundary, then 
Eq. (Q) represents an energy flow from the system. Between 
two times ti and ^2 one can integrate to find that Ail/ = 
— i Jg Sx^—^ T°-^%tlab is the change in the mass contained 
by Ctt- Subtraction terms from a reference spacetime do not 
need to be included here as it expresses the change in the mass 
of the system. The rate at which this work is done is 
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which describes the rate of change of the system's mass due 
to the purely gravitational interaction between it and the sur- 
rounding environment. 

It is illustrative to decompose the expression for the work 
into terms involving projections of $t7ah normal to and into 
the spatial two-surfaces Vlf. We find 
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where e = a kab/^T^, ja = cTabUcQ /Sir, and s 
[k"-'' + a°'''{n'^u'^AdUc — a'^^ kcd)] / S>7^ ai'e the quasilocal sur- 
face energy, momentum, and stress densities. The first two 
are potentials conjugate to changes in the lapse function and 
shift vector respectively while the surface stress density is a 
work potential conjugate to changes in the size and shape of 
the surface ilf 

The stress density can be further decomposed as follows. A 
change in the two-metric, Saab = QibSy/a + -Ja5<,ab, is writ- 
ten as a change in the "size" ^/a of the surface plus a change 
in the conformally-invariant part of the metric (the "shape" 
of the surface) = ^^abj \f^- Correspondingly, the surface 



stress density is decomposed into a surface tension s — s <^ab 
and a shear rf^ = s"^^ j \fa. Then we rewrite the work term as 
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The above has a particularly nice application in the physics 
of thin shells. Israel [|ll|] first showed that a thin shell of matter 
can be described in general relativity by matching two space- 
times along a timelike boundary B such that even though they 



induce the same surface metric on B, the extrinsic curvature 
in each spacetime is different. If 0^^^ and are those curva- 
tures this (mild) singularity can be accounted for if there is a 
(distributional) stress energy tensor Sab = '^ab^'''ab ^- ^ 
set of observers dwelling on the surfaces ilt (which foliate B) 
measures the shell to have matter-energy A/+ — [Eq. (^]. 
A more detailed discussion of this may be found in [||] but 
here we note that the above analysis for the quasilocal energy 
also shows that the set of observers dwelling on Vlt measures 
the matter-energy to change with rate dW/dt (Eq. Then 
the quasilocal densities defined above are the energy, angu- 
lar momentum, and stress tensor of the matter shell. A set 
of observers being evolved by = see work being done 
on the shell at a rate equal to the integral of the stress tensor 
contracted with the time rate of change of the area — exactly 
as one would expect from classical physics. 



m. GRAVITATIONAL RADIATION 

Equation (Q) purportedly measures the change in the mass 
of a system. In this section we apply our work formula to ob- 
tain the correct mass loss for a system radiating gravitational 
waves. For this we suppose that the quasilocal surface is in the 
wave-zone, far away from the radiating system. Although this 
is not a very interesting application of a quasilocal method 
(since an asymptotic method, such as the Bondi-Sachs mass 
loss formula could as well be used), it is useful to confirm that 
Eq. (Q does recover the correct result. 

Gravitational radiation far from the generating source 
can be described as a transverse-traceless perturbation to 
the flat-space metric. In spherical-polar coordinates, the 
metric is given by ds^ = —dt^ + dr^ + (rdd)^ + 
{r sined(l))^ + hf,^dx^dx'' where hf.^dx^^dx" = h+[{rdef - 
{r svaOd^iY'] + 2hx{rd9){r sin 9 dcf)) is the transverse, trace 
free perturbation. The "plus," h^, and "cross," h^, polar- 
izations represent outgoing, spherical waves, and have the 
form h+ {t, r, 6, (f)) — s+{t — r, 9, cf)) jr and /ix (t, f , ^, 0) = 
s-x{t ~ r,9,(f))/r. We then find the energy lost by the radiat- 
ing system by inserting this metric into Eq. (Q) while taking 
the boundary to be a sphere of constant r in the wave-zone 
(very large r). The integrand of Eq. (H) is 
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[idh+/dtf + idhx/dt)^ 
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to leading order in the perturbation and in r. This is the stan- 
dard expression for the flux of gravitational radiation — see, 
e.g.,Eq. (10) ofRef. 

By inspection of the form of the perturbation, it is clear 
that the energy loss arises due to the shearing of the bounding 
two-surface ilt since, to leading order, the perturbation does 
not affect the volume element on that two-surface. Thus the 
entire energy loss (in the transverse, trace-free gauge) arises 
from the "r]°'^$t'^ab' work term. 

As a simple example, consider two point-particles, each of 
mass m — M/ 2, orbiting each other in the xy-plane with an- 
gular frequency uj and constant separation a. The quadrupole 
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moment tensor X,7j in Cartesian coordinates is I^x = 
—lyy — iA/a'^ COS 2cjt and = iilfa^ sin 2tjt (constant 
terms omitted). The far-field metric perturbation is hjk = 
2(9^X,-fe/9t^)/r, so hyy = —h^x = (Afa^cj^/r) cos2cj(t — 
r) and h^y = — (A/a^w^/r) sin2[x)(i — r). 

Using Eqs. (4.3) and (4.4) of Ref. we find 

h+ = -^Ma^uo'^r-^ (I + cos^ 9) cos2[uo{t - r) - 4)] (7a) 
/ix = -A'/a^cj^r"^ cos6'sin2[a;(t- r) - 0]. (7b) 

We then integrate Eq. (^ over the sphere at large r to obtain 
the loss of energy from the system: 

- dE/dt^lM^a'^uu^ = l{M/af (8) 

where we have used Kepler's law a^iJ^ — M for particles in 
a circular orbit. 



IV. TIDAL HEATING 

We now calculate the work done by an external gravita- 
tional field to deform a self-gravitating body. The canonical 
example of this effect in the solar system is the tidal heating of 
lo by Jupiter. In this instance, the gradient of Jupiter's gravi- 
tational field distorts lo from being a perfect sphere and then 
tidally locks it in its orbit so that it always presents the same 
face to Jupiter That orbit is strongly perturbed by the other 
Gallilean moons and so its radial distance from Jupiter varies 
with time. With this variation comes a corresponding one in 
the gradient of the field and so lo is gradually stretched and 
then allowed to relax. The energy transferred by this pumping 
is largely dispersed as heat and it is this heat that produces the 
volcanic activity on lo. The same type of process occurs for 
any two bodies in non-circular orbits about each other 

First from a Newtonian perspective, we may mathemati- 
cally describe the gravitational fields in this situation as fol- 
lows. We assume that the self-gravitating body is far enough 
away from the source of the external field that that field is 
nearly uniform close to the body. Then in a rectangular co- 
ordinate system that orbits with the body with its origin at 
the center of mass, the Newtonian potential of the exter- 
nal field may be written as <i>ext = ^SijX^x^ where £ij is 
the (time-dependent but symmetric and trace-free) quadrupole 
moment of the field and a;* is the position vector based at the 
body's centre of mass. At the same time, to quadrupolar or- 
der the Newtonian potential of the body is $o = — A//r — 
^r^^Xijii'n^ , where A/ is the mass of the body, r is the radial 
distance from the centre of mass, is its (time-dependent 
but symmetric and trace-free) quadrupole moment, and ?i' is 
the unit normal radial vector. 

With this in mind the techniques of Thorne and Hartle 
can be used to construct a metric that describes these situa- 
tions in the slow moving, nearly Newtonian limit. First, de- 
fine an annulus surrounding the body whose inner boundary 
is chosen so that the gravitational field of the body is weak 
throughout and whose outer boundary is chosen so that the ex- 
ternal field is nearly uniform. This region is termed the buffer 



zone. The rectangular coordinate system is replaced with one 
that is chosen so that the metric is as close to Minkowskian as 
possible over the buffer zone [Q]. Then to first order in pertur- 
bations from Minkowski and first order in time derivatives the 
metric can be written as 

ds^ = -(1 + 2<P)dt^ + 2{Aj + dt^j)dx^dt 

+ [(1 - 2<^)5,j + + dj£.i\dx'dx' (9) 

where the indices run from one to three and 6ij is the Carte- 
sian metric diag[ 1,1,1] on a spacelike slice. The Newto- 
nian potential is $ = —M/r — ^{ir^^Iij — r^£ij)n^n^ and 
= -2r-'^n^dl,kldt - ^r^{bnjn^ - 25])n^d£ki/dt is 
a vector potential that must be added so that the metric is a 
solution to the first order Einstein equations. Here, n* is the 
radial normal with respect to the flat spatial metric 5ij and 

— x"^ + y"^ + z^. The diffeomorphism generating vector 
field represents the gauge ambiguity in setting up a nearly 
Minkowski coordinate system. In order that the metric be 
slowly evolving and nearly Minkowski, must be of the form 

— ar^'^Tjk'n'^ + Pr^Sjkn^ + 'jr^£kin'^n'"nj, where a, (3, 
and 7 are free constants of order one. 

We set up a constant r timelike quasilocal surface B in 
the buffer zone and foliate with constant t spacelike two- 
surfaces fit. Then the time vector is d/dt. In calculat- 
ing the rate of change of the mass contained within Q,t it 
is most convenient to switch to spherical coordinates. We 
make the standard transformation to spherical coordinates 

— r [sin 6 cos 0, sin sin (/), cos 0] ; in these coordinates, the 
metric is 

ds^ = -(1 + 2^)dt^ + (1 - 2$)[<ir2 + {rdef 
+ {rsined4>f] + 2Ardrdt + 2Ae{rde)dt 
+2A^{rsin0d(j))dt + Hrrdr'^ + Hee{rdef 
+H^^{rs\ned(t)f + 2Hredr{rde) 
+2Hr^dr{r sin ed<j)) + 2Hg^{rd6){r sin ed<j)f (10) 

where Hrr — —Aar^^Irr + 6(/? + "f)r'^£rr, Hgg — 

2ar-^Iee + 2(3r'^£ge + 2-ir^£rr, H^^ = 2ar-^I^^ + 
2(3r^£^^ + 2-fr'^£rr, H^e = -ar-^Ire + (4/3 + 2jy£re, 
Hr^ = -ar-^Ir^+{AP+2j)r^£r4„ and Hg^ = 2ar-^Ie^+ 
2(3r'^£g^. In these expressions £„ = £ije\.el, £rg = £ije].eg, 
etc., with e* = n', e\ — dge\. and — [1/ sin9)d^e\.. Also, 

Ar ~ [Aj + dtij)el, etc., but we don't need their expanded 
forms since only time derivatives of them show up in later cal- 
culations and we are ignoring second order time derivatives. 

As might be expected, the subsequent calculations are quite 
involved and we did them partially with GRTensor [|l4[]. To 
lowest order 

dW I f ,2 , „h^ 

+47^ + 8/37)r5£,,5„- 

+2(3 - 2a + 6/3 - I27 + 8aj)£,jl^j 

-(-9 + 12a + 4Q2)r-%X,j]. (11) 
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The calculations used the identities J^^ d9d(j) 
{8tt/ 15) AijBij and J^^^ d6d(l)sin0{2Ae4,Be4, - AeeB^^ - 
AijjijjBoe) — (47r/3)Ayi?y where the integrations are over 
the unit sphere. 

This result requires some interpretation. As the external 
field changes with time and thereby forces the self-gravitating 
body to change configuration, the work done by the external 
field can be split into time reversible and irreversible parts [as 
seen in Eq. (|ll|)]. The reversible work represents work be- 
ing done to increase the potential energy of the system and 
is recoverable. On the other hand the irreversible part repre- 
sents work done to deform and/or heat up the system. This 
is the tidal heating that we are interested in. Further, from the 
quasilocal perspective, we expect to see an energy flow arising 
from fluctuations of the quasilocal surface within otherwise 
static fields. Of course this work would also be reversible. 
Thus, it is only the irreversible part that we are interested in 
and we have calculated that to be ^Eijdlij / dt above. This 
is the same leading term obtained when one does the corre- 
sponding calculation in Newtonian gravity or with pseudoten- 
sors and it is independent of diffeomorphisms generated 
by which correspond to fluctuations of the quasilocal sur- 
face. Note however, that the time reversible and gauge de- 
pendent terms of equation ( 11 ) are dependent on those fluctu- 
ations and furthermore that dependence is different from that 
found in ref. |Q| using pseudo-tensor methods. Similarly other 
pseudo-tensor or quasilocal methods would obtain a different 
gauge dependence for these terms. What is important is that 
the physically relevant time irreversible term does not depend 
on the -generated diffeomorphisms. 

Finally for completeness let us consider how this en- 
ergy flow splits up into its components parts as consid- 
ered in Eq. (Q). Then to the order that we are interested 
the angular momentum term is zero and we are left with 
two terms dWn/dt = - J d9d(f)y/ae$tN and dWa/dt = 
\ J ded(j)^Ns''^%tcjab. We find 
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dt 2 dt 



a dSij _ /3 dirj 
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15 dt 



dt 
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-3(2a-3)rt„Xy]. 
The second term is a bit more complicated. It is 

dWrr a dSi i „ 6 „ dl, 



dt 



(12) 



l^iX - 

dt 15 dt 

1 d 



0. 
5' 



- dt 



47 dljj 
5 dt 



+ ^^[(-1 - 3/3 - 2/3^ + 47^ + SP^yS^jE,, 
+2(3 - a + 3/3 - 67 + Aa-/)£,jJij 
-(2a'-9a + 9)r-%2'y]. (13) 

Thus part of the work done is measured by deformations of 
the surface and part is measured by changes in how observers 
choose to measure the rate of passage of time. Note that in- 
dividually the time irreversible sections of the two parts are 
gauge dependent, but when we combine them we reobtain 



Eq. (11) and the gauge dependence vanishes back into the re- 
versible part where we would expect it. 



V. CONCLUSIONS 

We have modified the quasilocal energy formalism of 
Brown and York so that it may be used to study non-stationary 
spacetimes where energy flows in and out through the quasilo- 
cal surface. As applications of this extension we have exam- 
ined implications for the physics of relativistic thin shells of 
matter, the energy carried from a source to infinity by grav- 
itational waves, and the transfer of energy to a body during 
gravitational tidal heating. The success of the formalism in all 
three applications provides further evidence that the Brown- 
York energy has physical content. Furthermore, in the tidal 
heating application we have seen how the quasilocal formal- 
ism provides a geometrical explanation of the gauge ambigu- 
ities that are also found in the Newtonian and pseudotensor 
approaches. 
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